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ABSTRACT 



A review on acoustic mode damping and excitation in 
solar-type stars is presented. Current models for lin- 
ear damping rates are discussed in the light of recent 
low-degree solar linewidth measurements with emphasis 
on the frequency-dependence of damping rates of low- 
order modes. Recent developments in stochastic exci- 
tation models are reviewed and tested against the latest 
high-quality data of solar-like oscillations, such as from 
alpha Cen A, and against results obtained from hydrody- 
namical simulations. 

Key words: mode physics; mode damping, stochastic ex- 
citation. 



1. INTRODUCTION 

Solar p modes are believed to be excited stochastically in 
the outer layers of the star by the turbulent convection. 
Accurate measurements are available from both space- 
borne (e.g., VIRGO, GOLF and MDI) and ground-based 
(GONG and BiSON) observations which provide valu- 
able data regarding the excitation and damping processes 
influencing the modes. Such modes can be stochasti- 
cally excited by the turbulent convection. The process 
can be regarded as multipole acoustical radiation (e.g. 
Unno 1964). For solar-like stars, the acoustic noise gen- 
erated by convection in the star's resonant cavity may be 
manifest as an ensemble of p modes over a wide band 
in frequency (Goldreich & Keeley 1977). The ampli- 
tudes are determined by the balance between the exci- 
tation and damping, and are expected to be rather low. 
The turbulent-excitation model predicts not only the right 
order of magnitude for the p-mode amplitudes (Gough 
1980), but it also explains the observation that millions of 
modes are excited simultaneously. One particular obser- 
vational detail, the frequency dependence of the energy 
supply rate to acoustic modes, has been shown to be a 
particularly important diagnostic property (e.g. Libbrecht 
1988), for it has diagnostic potential for improving cur- 
rent stochastic-excitation models (e.g., Balmforth 1992b, 
ChapHn et al. 2005). 

The energy flow from radiation and convection into and 
out of the p modes takes place very near the surface. Sun- 
like stars possess surface convection zones, and it is in 



these zones, where the energy is transported principally 
by the turbulence, that most of the driving takes place. 
Mode stability is governed not only by the perturbations 
in the radiative fluxes (via the K-mechanism) but also by 
the perturbations in the turbulent fluxes (heat and momen- 
tum). The study of mode stability therefore demands a 
theory for convection that includes the interaction of the 
turbulent velocity field with the pulsation. 

Our understanding of turbulent convection in stars is still 
only rudimentary. Simple phenomenological convection 
models based on mixing-length ideas are still the most 
widely used for computing just the mean stratification of 
the convectively unstable layers. Attempts to describe 
the associated physics on all scales (both temporal and 
spatial) have led to similar simplified formalisms, which 
include one or more adjustable parameters that require 
calibration by high-quality observations. 

2. MODE PARAMETERS 

Were solar p modes to be genuinely linear and sta- 
ble, their power spectrum could be described in terms 
of an ensemble of intrinsically damped, stochastically 
driven, simple-harmonic oscillators, provided that the 
background equilibrium state of the star were indepen- 
dent of time (Fig. 1); if we assume further that mode 
phase fluctuations contribute negligibly to the width of 
the spectral lines, the intrinsic damping rates of the 
modes, r/2, could then be determined observationally 
from measurements of the pulsation linewidths F. 

The power (spectral density), P, of the surface displace- 
ment (,ni{t) of a damped, stochastically driven, simple- 
harmonic oscillator, satisfying 
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which represents the pulsation mode of order n and de- 
gree /, with linewidth r„; and frequency ujni and mode 
inertia I^i, satisfies 
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assuming Tni <C uJni, where f{t) describes the stochastic 
forcing function. Integrating equation dU over frequency 
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Figure 1. Power spectral density of a randomly excited, 
damped, harmonic oscillator Pf represents the spectral 
density of the random force and P is the product of the 
Lorentzian and Pf (adopted from Kosovichev 1995). 

leads to the total mean energy in the mode 



P{uj) duj 
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(3) 



where Ani is the displacement amplitude (angular brack- 
ets, (), denote an expectation value) and Vnj is the rms 
velocity of the displacement. The total mean energy of a 
mode is therefore directly proportional to the rate of work 
(also called excitation rate or energy supply rate) of the 
stochastic forcing Pf at the frequency ujni and indirectly 
proportional to Tni- 

Equations (dJ-© are discussed in terms of the displace- 
ment ^ (from here on we omit the subscripts n and I), but 
in order to have a direct relation between the observed 
velocity signal v{t) = d£_/dt and the modelled excitation 
rate Pf we shall first take the Fourier transform ¥{1/) of 
v{t) {v = uj/2n). It follows that the total mean energy E 
of the harmonic signal of a single pulsation mode is then 
given by (Chaplin et al. 2005) 



E = IV= Id 



\Viiy)fdiy=^ITH, (4) 



(5) 



is the maximum power density - which corresponds to 
the 'height' of the resonant peak in the frequency domain 
(see Fig. 1). The height H is the maximum of the dis- 
crete power, i.e. the integral of power spectral density 
over a frequency bin 5 ~ l/Tobs, where Tobs is the total 
observing time. The following expressions 



and 
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H := 



(r/2)2/ 
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(7) 



provide a direct relation between the observed height H 
(in cm^ s~^Hz~^), the modelled energy supply rate Pf (in 
erg s~^), and damping rate 77 = r/2. 



3. MODEL COMPUTATIONS 

The model calculations need to solve the fully nonadia- 
batic pulsation problem from which the damping rates are 
obtained from the imaginary part r/ = of the complex 
eigenfrequency cj = tJr+it^i- Moreover, the pulsation dy- 
namics depends crucially on the convection treatment as 
reported by Balmforth (1992a) for the solar case and by 
Houdek et al. (1999) for solar-type stars. Consequently 
the model computations should also take into account the 
effect of the turbulent fluxes (heat and momentum) on the 
pulsation properties (e.g., damping rates and shape of the 
eigenfunctions). 

Here we follow the calculations reported by Balmforth 
(1992a) and by Houdek et al. (1999). In these calcu- 
lations the equations describing the stellar structure and 
pulsations are 
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resembling the conservation of momentum, mass and 
thermal energy. In equations (O is r the radius, m is 
mass, pg is gas pressure, p is density, T is tempera- 
ture, Cp is the specific heat at constant pressure, and 
S := -~{d\n p/d\nT)p . The turbulent velocity field 
u = (mi,M2,M3) is described in a cartesian coordinate 
system and for a Boussinesq fluid, pt = {pu'^) is the 
(r, r) -component of the Reynolds stress tensor (some- 
times called the turbulent pressure, and here the angular 
brackets denote an ensemble average). In the Boussi- 
nesq approximation the convective heat flux is ~ 
pcp{u3-d) being the convective temperature fluctua- 
tion), and ^ = {u ■ u) / (uf) is a parameter that describes 



the anisotropy of the turbulent velocity field u. The ra- 
diative heat flux F,- is obtained from the Eddington ap- 
proximation to radiative transfer (Unno & Spiegel 1966). 
The equilibrium model is constructed from solving equa- 
tions (|8]l but with the partial time derivatives set to zero, 
d/dt = 0, which, for example, leads for the first integral 
of the thermal energy equation to the well-known expres- 
sion L{r) = 47rr2(Fi. + Fc). 

The pulsation equations for radial modes are then ob- 
tained from linear perturbation theory 
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where S is the Lagrangian perturbation operator, and for 
simplicity the right hand side of the perturbed momen- 
tum equation is formally expressed by the function / (the 
full set of equations can be found in, e.g., Balmforth 
1992a). Equations (|9]l are solved subject to boundary 
conditions to obtain the eigenfunctions and the complex 
angular eigenfrequency ui — LUy + irj, where uj,- is the 
(real) pulsation frequency and rj — r/2 is the damping 
rate in (s^^). The turbulent flux perturbations of heat and 
momentum, SLc and Spt, and the fluctuating anisotropy 
factor (5$ are obtained from the nonlocal, time-dependent 
convection formulation by Gough (1977a,b). 



4. AMPLITUDE RATIOS AND PHASES IN THE 
SOLAR ATMOSPHERE 

A useful test of the pulsation theory, independent of an 
excitation model, can be performed by comparing the the- 
oretical intensity-velocity amplitude ratios 



AV 



SL/L 
uji-r 5r/r 



(10) 



with observations. It is, however, important to realize 
that various instruments observe in different absorption 
lines and consequently at different heights in the atmo- 
sphere. This property has to be taken into account not 
only when comparing observations between various in- 
struments (e.g. Christensen-Dalsgaard & Gough 1982, 
see also Bedding, these proceedings), but also when com- 
paring theoretical amplitude estimates with observations 
(Houdek et al. 1995). Table 1 lists some of the relevant 
properties of various instruments. 

In the top panel of Fig. |2]the theoretical amplitude ratios 
(equation ( fTOl l) of a solar model are plotted as function of 
height for several radial pulsation modes. The mode en- 
ergy density (which is proportional to rp^^^Sr) increases 
rather slowly with height; the density p, however, de- 
creases very rapidly and consequently the displacement 
eigenfunction Sr increases with height. This leads to the 



Table 1. Absorption lines and their wavelengths A of var- 
ious helioseismic instruments. Also listed are the optical 
depths T5000 at 5000 A and the corresponding approxi- 
mate heights above the photosphere (h — at T — T^s) 
at which the lines are formed. 



Instrument 


line 


A (A) 


T'SOOO 


height (km) 


BBSO 


Ca 


6439 


0.05 


- 129° 


BiSON 


K 


7699 


0.013 


- 250'' 


MDI 


Nil 


6708 


9x10- 


3 - 300= 


GOLF 


NaDl/D2 


5690 


5x10- 


- 500= 



'^fromLibbrecht(1988), 
^from Christen-Dalsgaard & Gough (1982), 
=from Toutain et al. (1997, but see also 
Baudin et al., these proceedings). 

results shown in the upper panel of Fig. |2] where the de- 
crease in the amplitude ratios with height is particularly 
pronounced for high-order modes for which the eigen- 
functions vary rapidly in the evanescent outer layers of 
the atmosphere. It is for that reason why solar velocity 
amplitudes from, e.g., the GOLF instrument have larger 
values than the measurements from the BiSON instru- 
ment (by about 25%, Kjeldsen et al 2005). 
The lower panel of Fig. |2] compares the estimated so- 
lar amplitude ratios (curves) with observed ratios (sym- 
bols) as function of frequency. The model results are de- 
picted for velocity amplitudes computed at different at- 
mospheric levels. The observations are obtained from 
accurate irradiance measurements from the IPHIR in- 
strument of the PHOBOS 2 spacecraft with contempo- 
raneous low-degree velocity data from the BiSON instru- 
ment at Tenerife (Schrijver et al. 1991). The thick solid 
curve represents a running-mean average, with a width of 
300 /iHz, of the observational data. The theoretical ratios 
for h = 200 km (dashed curve) show reasonable agree- 
ment with the observations. 

An additional useful property of the oscillations is 
the phase difference between intensity and velocity. 
Phase shifts between intensity observations from the 
VIRGO/SPM instrument and contemporaneous velocity 
measurements from the GOLF instrument (symbols with 
error bars) are plotted as function of frequency in Fig. [5] 
together with results from various model computations 
(curves). The dashed curves are the results from cal- 
culations by Marmolino & Severino (1991) who treated 
the energy transport with Newton's cooling law assum- 
ing various values for the damping time; the dot-dashed 
curve is the result by Gough (1985) who used the diffu- 
sion approximation to radiative transfer and a local, time- 
dependent mixing-length model to estimate the turbulent 
fluxes. The solid curves are the results from the model 
discussed in Section 3 and for three different sets of con- 
vection parameters (Houdek et al. 1999). The latter re- 
sults are in good agreement with the observations. 
The results depicted in Figs |2] and [3] provide us with 
confidence that the pulsation computations reproduce the 
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Figure 3. Intensity-velocity phase differences of low- 
degree p modes at the maximum of solar activity. Symbols 
are the results from intensity observations in the green 
channel of the VIRGO/SPM instrument on the SOHO 
spacecraft, and from contemporaneous velocity obser- 
vations from the GOLF instrument. Theoretical results 
are depicted for three models: the dashed curves assume 
Newton's cooling law for energy transfer with three dif- 
ferent values for the damping times, the dot-dashed curve 
includes a local treatment of convection dynamics and 
the diffusion approximation to radiative transfer, and the 
solid curves are the results from the model discussed in 
Section 3 for three different sets of convection parame- 
ters (adopted from Jimenez 2002). 



Figure 2. Top: Calculated amplitude ratios ( see equa- 
tion AlOi j as a function of height in a solar model for 
modes with different frequency values. Bottom: Theo- 
retical amplitude ratios ( surface luminosity perturbation 
over velocity) for a solar model compared with observa- 
tions by Schrijver etal. ( 1991 ). Computed results are de- 
picted at different heights above the photosphere (h=0 km 
at T — Tcs) The thick, solid curve indicates a running- 
mean average of the data (from Houdek et al. 1995). 

properties of the radial eigenfunctions in the outer atmo- 
spheric solar layers reasonably well. 

5. SOLAR DAMPING RATES 

Damping of stellar oscillations arises basically from two 
sources: processes influencing the momentum balance, 
and processes influencing the thermal energy equation. 
Each of these contributions can be divided further accord- 
ing to their physical origin, which was discussed in detail 
by Houdek et al. (1999). 

Important processes that influence the thermal energy 
balance are nonadiabatic processes attributed to the mod- 
ulation of the convective heat flux by the pulsation. This 
contribution is related to the way that convection mod- 
ulates large-scale temperature perturbations induced by 
the pulsations which, together with the conventional k- 
mechanism, influences pulsational stability. 
Current models suggest that an important contribution 
that influences the momentum balance is the exchange 



of energy between the pulsation and the turbulent ve- 
locity field through dynamical effects of the fluctuating 
Reynolds stress. In fact, it is the modulation of the turbu- 
lent fluxes by the pulsations that seems to be the predom- 
inant mechanism responsible for the driving and damp- 
ing of solar- type acoustic modes. It was first reported 
by Gough (1980) that the dynamical effects arising from 
the turbulent momentum flux perturbation dpt contribute 
significantly to the damping F. Detailed analyses (Balm- 
forth 1992a) reveal how damping is controlled largely by 
the phase difference between the momentum perturbation 
and the density perturbation. Therefore, turbulent pres- 
sure fluctuations must not be neglected in stability analy- 
ses of solar-type p modes. 

A comparison between the latest linewidth measurements 
(full- width at half-maximum) A„; — Tni/^TT and the- 
oretical damping rates is given in Fig. |4] The observa- 
tional time series from BiSON (Chaplin et al. 2005) was 
obtained from a 3456-d data set and the linewidths of the 
temporal power spectrum (top panel) extend over many 
frequency bins 6 = l/2robs- In that case the linewidth 
in units of cyclic frequency is related to the damping rate 
according to 

A„( ^ 7r^i?7„i . (11) 

Recently Dupret et al. (2004, see also these proceed- 
ings) performed similar stability computations for the 
Sun using the time-dependent mixing-length formulation 
by Gabriel et al. (1975, 1998) which is based on the for- 
mulation by Unno (1967). The outcome of their com- 
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Figure 4. Top: Power spectrum of solar low-degree p modes obtained from a 3456-d data set collected by BiSON between 
1991 and 2000 ( Chaplin et al. 2005). The two insets show Lorentzian profile fits (and their full-width at half-maximum 
(FWHM); solid, blue curves) to the spectral peaks of radial modes with order n = 13 (left) and n = 21 (right); both 
spectral peaks are indicated by vertical arrows in the power spectrum. Bottom: The symbols in the left-hand panel 
are the measured BiSON linewidths A„; = Tni/^n (we denote the FWHM in unit of cyclic frequency by A„;j which are 
compared with the theoretical damping rates tt^ ^ rjni ( connected by the solid curve ) obtained from the model computations 
of Section 3 (from Chaplin et al. 2005). In the right-hand panel theoretical results of {2Tr)^^r]ni (solid curve) by Dupret 
et al. (2004) are compared with observations o/A„;/2 (symbols). 



putations is illustrated in the lower right panel of Fig. |4] 
which also shows the characteristic plateau near 2.8 mHz. 
It is, however, interesting to note that their findings sug- 
gest the fluctuating convective heat flux to be the main 
contribution to mode damping, whereas for the model re- 
sults shown in the lower left panel of Fig. 2] which are 
based on Gough's (1977a,b) convection formulation, it is 
predominantly the fluctuating Reynolds stress that makes 
all modes stable. 



6. EXCITATION MODEL AND SOLAR AMPLI- 
TUDES 

Because of the lack of a complete model for convection, 
the mixing-length formalism still represents the main 
method for computing the turbulent fluxes in the convec- 
tively unstable layers in a star. One of the assumptions in 
the mixing-length formulation is the Boussinesq approx- 
imation, which results in neglecting the acoustic wave 



generation by assuming the fluid to be incompressible. 
Consequently a separate model is needed to estimate the 
rate of the acoustic noise (energy supply rate) generated 
by the turbulence. The excitation process can be regarded 
as multipole acoustic radiation (Lighthill 1952). Acous- 
tic radiation by turbulent multipole sources in the con- 
text of stellar aerodynamics has been considered by Unno 
& Kato (1962), Moore & Spiegel (1964), Unno (1964), 
Stein (1967), Goldreich & Keeley (1977), Bohn (1984), 
Osaki (1990), Goldreich & Kumar (1990), Balmforth 
(1992b), Goldreich, Murray & Kumar (1994), Musielak 
et al. (1994), Samadi & Goupil (2001) and Chaplin et al. 
(2005). 

6.1. Excitation model 

The mean amplitude A of a mode is determined by a bal- 
ance by the energy supply rated Pf from the turbulent 
velocity field and the thermal and mechanical dissipation 
rate characterized by the damping coefficient rj. The pro- 



cedure that we adopt to estimate A is that of ChapHn 
et al. (2005), whose prescription follows that of Balm- 
forth (1992b). 

As in Section 2 we represent the linearized pulsation dy- 
namics by the simplified equation 
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g + 2,| + /:|)=^N + g(,s') 



(12) 



for the displacement ^(r , t), which is now also a function 
of radius r, of a forced oscillation corresponding to a sin- 
gle radial mode satisfying the homogeneous equation 



(13) 



in which uj (and ^) are real and £ is a linear spatial op- 
erator. The (inhomogeneous) fluctuating terms on the 
right-hand-side of equation (fT2l) arise from the fluctuat- 
ing Reynolds stresses 



^{u) = V • {puu — {puu)) 



(14) 



and from the fluctuating gas pressure (due to the fluctu- 
ating buoyancy force), represented by Q{s'), where s' is 
the Eulerian entropy fluctuation (Bohn 1984; Osaki 1990; 
Goldreich & Kumar 1990; Balmforth 1992b; Goldreich, 
Mun-ay & Kumar 1994; Samadi & Goupil 2001). The 
latest numerical simulations by Stein et al. (2004) sug- 
gest that both forcing terms in equation (fTSI) contribute 
to the energy supply rate Pf by about the same amount, 
a result that was also reported by Samadi et al. (2003) 
using the turbulent velocity field and anisotropy factors 
from numerical simulations (Stein & Nordlund 2001). In 
this paper we consider only the term of the fluctuating 
Reynolds stresses and because we use only radial modes, 
only the vertical component of is important. 



d_ 

dr 



{puD). 



(15) 



If we define the vertical component of the velocity cor- 
relation as i?33 its Fourier transform R^s can 
be expressed in the Boussinesq-quasi-normal approxima- 
tion (e.g. Batchelor 1953) as a function of the turbulent 
energy spectrum function E{k, oj): 



R. 



33 



127rfc2 



(16) 



where fc is a wavenumber and 5* is an anisotropy param- 
eters given by 



* = [2$/3(4>- 1)] 



1/2 



(17) 



which is unity for isotropic turbulence (Chaplin et al. 
2005). This factor was neglected in previously published 
excitation models but it has to be included in a consistent 
computation of the acoustic energy supply rate. Follow- 
ing Stein (1967) we factorize the energy spectrum func- 
tion into E{k,uj) = E{k)n{uj;Tk), where = X/kuk 
is the correlation time-scale of eddies of size n/k and ve- 
locity lifc; the correlation factor A is of order unity and 
accounts for uncertainties in defining r^. 
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Figure 5. Comparison of the frequency factor il,{uj;Tk) 
(see equation M9h between hydrodynamical simulation 
results for solar convection (Stein & Nordlund 2001; 
symbols) and the Lorentzian frequency factor by Samadi 
et al. (2003) for two different values of the correlation 
factor A (solid and dashed curves). The results are plot- 
ted as function of frequency evaluated 0.99 Mm below the 
solar photosphere (adopted from Samadi et al. 2003). 



The energy supply rate is then given by (see Chaplin et al. 
2005 for details) 



Pf 



91 



^3 I ^q,j.p^ 



S(r;Lu)dr, (18) 



with 



K-^E^{K)n{Tk;Lu)dK, (19) 



where k = ki/n, £ is the mixing length, R is surface ra- 
dius, and is the normalized radial part of ^. The spec- 
tral function S accounts for contributions to Pf from the 
small-scale turbulence and includes the normalized spa- 
tial turbulent energy spectrum E{k) and the frequency- 
dependent factor fl{Tk;uj). For E{k) it has been com- 
mon to adopt either the Kolmogorov (Kolmogorov 1941) 
or the Spiegel spectrum (Spiegel 1962). The frequency- 
dependent factor fl{Tk ; oj) is still modelled in a very rudi- 
mentary way and we adopt two forms: 
- the Gaussian factor (Stein 1967) , 



(20) 



- the Lorentzian factor (Gough 1977b; Samadi et al. 
2003; Chaplin et al. 2005), 



1 



7r\/2 In 2 1 + (cjTfc/\/2 In 2)2 



(21) 



The Lorentzian frequency factor is a result predicted for 
the largest, most-energetic eddies by the time-dependent 
mixing-length formulation of Gough (1977b). Recently, 
Samadi et al. (2003) reported that Stein & Nordlund's hy- 
drodynamical simulations also suggest a Lorentzian fre- 
quency factor (see Fig. |5]l, which decays more slowly 




Depth (Mm) 



Figure 6. Top: The contours represent the integrand 
of the acoustic energy supply rate (\18h multiplied by 
I^^dr /dlnp, and are plotted as function of frequency 
V and the depth variable z = Rq— r. The plus symbols 
with ( horizontal) error bars are the measured locations of 
the excitation regions by Chaplin & Appourchaux ( 1999). 
The solid and dashed curves denote the exact acoustic 
cutoff frequency and the approximation in an isothermal 
region, cjlHp, where c is the sound speed and Hp is the 
pressure scale height. Bottom: Logarithm of integrand of 
the acoustic energy supply rate from Stein & Nordlund's 
(2001 ) hydrodynamical simulations. 

with depth z and frequency uj than the Gaussian factor. 
Consequently a substantial fraction to the integrand of 
equation (fTSl l arises from eddies situated in the deeper 
layers of the Sun, resulting in a larger acoustic excitation 
rate Pf. 

Another way to study the properties of the turbulent- 
excitation model is to compare the depth and frequency 
dependence of the integrand of Pi with solar measure- 
ments and with hydrodynamical simulations of the Sun. 
The results are depicted in Fig.|6] the contours in the top 
panel show the integrand of equation ( fTSl l multiplied by 
/^^dr/dlnp. The contours are compared with measure- 
ments of the locations of the driving regions by Chaplin 
& Appourchaux's (1999), represented by the plus sym- 
bols with (horizontal) error bars. The modelled loca- 
tions of the excitation regions are in good agreement with 
the observations, particularly for modes with frequencies 
larger than the acoustic cutoff frequency (these modes are 
therefore propagating). The cutoff frequency is the solid 
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Figure 7. Reynolds stress as a function of the depth vari- 
able z = Rq — r for various solar models. Results 
are shown for the nonlocal mixing-length model ( solid, 
black curve) and from hydrodynamical simulations by 
Trampedach et al. (1999, dashed, red curve) and Lud- 
wig (2005, dot-dashed, blue curve). 



curve and the dashed curve is the isothermal approxima- 
tion c/2Hp, where c is the sound speed and Hp is the 
pressure scale height. For these (propagating) modes the 
frequency dependence of Pi is predominantly determined 
by the turbulent spectrum (represented by the spectral 
function S, cf. equation (fT9ll). Modes with frequencies 
smaller than the acoustic cutoff frequency are evanescent, 
and their frequency dependence of Pi is predominantly 
determined by the shape of the eigenfunctions, i.e. by the 
term {d^r / drY / 1 , and consequently by the structure of 
the equilibrium model. 

The lower panel of Fig.|6]shows the logarithm to the base 
10 of the integrand of the excitation rate of the numeri- 
cal simulations by Stein & Nordlund (2001). As for the 
model results (top panel), the extent of the driving regions 
decreases with frequency; however, the simulations show 
broader excitation regions than the model results. This is 
predominantly a result of the different spatial extents of 
the Reynolds stresses {pu^ between the hydrodynami- 
cal simulations and the nonlocal mixing-length model. A 
plot of (/3u§) Ip [p is the total pressure) as a function of 
the depth variable z — Kq — r is illustrated in Fig. |7] 
The Reynolds stress of the nonlocal mixing-length model 
shows a narrow peak near z ~ 120 km and falls off more 
rapidly with z than the results from both hydrodynami- 
cal simulations. This leads to a smaller excitation rate 
for the mixing-length model compared with the hydrody- 
namical results and consequently the modelled heights H 
need to be scaled with a scaling factor A > 1 in order to 
reproduce the observed values of the mode peak heights 
(Chaplin et al. 2005). 

6.2. Solar amplitudes 

Model predictions of the mode height H were computed 
according to equation ^ assuming $ to be constant; the 
value 1.13 was adopted because it produces the best fit 
to the linewidth data (cf. lower left panel of Fig. |4]i. 
The resulting values are joined by lines of various style 




1000 



sooo 




I 3032 30C 

iilljiliM, 



3000 

Frequency (/J-Hz) 



4000 



5000 



10^ 



7 10" 



10= 



Kolmogorov: X Spiegel; \ 

0.3 0.3 

1.0 1.0 

1.7 1.7 




T i i + BiSON data 
X model results scaled to max(data) 



1000 



lO^i 



1500 2000 2500 3000 3500 4000 
I' (mHz) 



Kolmogorov: X Spiegel: A 

3 0.3 

1.0 1.0 

1.7 1.7 




1000 1500 2000 2500 3000 3500 4000 



Figure 8. Top: Power spectrum of solar low-degree p modes obtained from a 3456-d data set by BiSON (Chaplin et 
al. 2005). The two insets show Lorentzian profile fits (and their heights H; solid, blue curves) to the spectral peaks of 
radial modes with order n = 13 (left) and n = 21 (right); both spectral peaks are indicated by vertical arrows in the 
power spectrum. Bottom: Radial mode peak heights, H, from observation and modelling. The crosses (with error bars) 
are well determined heights extracted from fits to the 3456-d BiSON spectrum. The curves show predicted heights H for 
two turbulent spectra, due to Kolmogorov and Spiegel, with a Gaussian factor Q.q, and assuming various values for the 
eddy correlation factor A; they are scaled such that their maximum values agree with the maximum value of the data. 
Left-hand panel - the theoretical curves were computed using modelled values for both the forcing Pf and linear damping 
rates rj. Right-hand panel - the theoretical curves were computed using the modelled values for the forcing Pf but with 
the damping rates replaced by tt A, where A are observed mode widths. 



in the lower panels of Fig. [8] for different values of the 
correlation factor A and for the two different forms of 
the spatial energy spectrum E{k), in all cases with the 
Gaussian frequency factor given by equation ( l20b . The 
BiSON measurements are rendered as crosses (with er- 
ror bars). The maximum values of the theoretical results 
have been scaled to the maximum observed value. The 
values of the scaling factor A by which we have mul- 
tiplied the theoretical values are given in Chaplin et al. 
(2005). The scaled theoretical estimations of H shown in 
the lower left-hand panel of the figure used modelled val- 
ues of both the energy supply rate Pf and the linear damp- 
ing rate ?/. They exceed the observed values at both low 
and high frequency, largely as a consequence of the mod- 
elled linewidths being larger than their measured coun- 
terparts at intermediate frequencies. This is demonstrated 
by comparison with the lower right-hand panel of Fig. [8] 
where predictions of H have been made using the mod- 



elled energy supply rate but with the observed BiSON 
linewidths A, rather than the modelled damping rates rj. 
There is a dramatic improvement in the modelled spec- 
tral heights H, although agreement with observations re- 
mains imperfect. The discrepancy that remains when us- 
ing observed linewidths is indicative of error in the equi- 
librium model and in the modelling of the energy supply 
rate Pf. This comparison reveals that the principal error 
in the theory, at least for the Sun, rests in the calculation 
of the damping rates, particularly for low-order modes. 
These modes have upper turning points well below the 
region where most of the excitation takes place. 
Chaplin et al. (2005) have also computed mode heights 
H using the Lorentzian frequency factor (|2TI) . but have 
found that in that case the heights at low frequency are 
severely overestimated. This result comes about be- 
cause the Lorentzian factor decreases too slowly with 
depth at constant frequency. This is particularly so for 
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Figure 9. Top: Predicted apparent velocity amplitudes 
(defined to be times the rms value) for a model of 
a Cen A, computed according to equation Bottom: 
Linear damping rates for a model of a Cen A, obtained 
by solving the fully nonadiabatic pulsation equations dP]). 

modes of lowest frequency, for which the Lorentzian is 
widest. Since it is not obvious where in the frequency 
spectrum the transition from a Lorentzian (largest most- 
energetic eddies) to a Gaussian (small-scale turbulence) 
frequency spectrum takes place, results are shown only 
for the purely Gaussian frequency factor 

7. SOLAR-LIKE AMPLITUDES IN OTHER 
STARS 

Fairly accurate measurements of solar-like oscillation 
amplitudes in other stars are available today from ground 
based observations (e.g. Kjeldsen et al. 2005; Bed- 
ding, these proceedings). Estimates of solar-like os- 
cillation amplitudes in other stars were carried out by 
Christensen-Dalsgaard & Frandsen (1983), Kjeldsen & 
Bedding (1995), Houdek et al. (1999), Houdek & Gough 
(2002) and Samadi et al. (2005). 

As in the solar case, damping rates for other stars are ob- 
tained from solving the eigenvalue problem (|9]l and from 
calculating the excitation rates from expression ( fTSI ). 
With these estimates for -q and Pf Houdek & Gough 
(2002) predicted velocity amplitudes for several stars, us- 
ing equation (|6]l. Results for stochastically excited oscil- 
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Figure 10. Top: Predicted apparent velocity amplitudes 
(defined to be \/2 times the rms value) for a model of 
£, Hydrae, computed according to equation dS)). Bottom: 
Linear damping rates for a model of ^ Hydrae, obtained 
by solving the fully nonadiabatic pulsation equations (O 
(adopted from Houdek & Gough 2002). 



lation amplitudes and linear damping rates in the solar- 
like star a Cen A and in the sub-giant ^Hydrae are il- 
lustrated in Figs|9]and[T0l Bedding etal. (2004) reported 
mode lifetimes for a Cen A between 1-2 days (see 
also Fletcher, these proceedings) which are in reason- 
able agreement with the theoretical estimates of about 1.7 
days for the most prominent modes (the mode lifetime 
T — rj^^; see lower panel of Fig. |9l). For ^ Hydrae, how- 
ever, the theoretical mode lifetimes of the most prominent 
modes are t ~ 17 days which are in stark contrast to the 
measured values of about 2-3 days by Stello et al. (2004, 
2006), yet the estimated velocity amplitudes for ^ Hydrae 
are in almost perfect agreement with the observations by 
Frandsen et al. (2002). The estimated velocity ampli- 
tudes, together with theoretical results for other solar-like 
oscillators, are compared with measurements from vari- 
ous observing campaigns and with the suggested scaling 
laws by Kjeldsen & Bedding (1995, 2001) in Table 2. 
For the cooler stars the theoretical results are in reason- 
able agreement with the observations. For the rather hot- 
ter star Procyon, however, the theoretical velocity ampli- 
tudes are severely overestimated (see Table 2). 
This comparison between predicted and observed veloc- 
ity amplitudes is illustrated in Fig. [TT] The dotted line 



Table 2. Comparison of scaling laws and theoretical velocity estimates with the observations of solar-like oscillations. 
The quantity V^bs is the typical observed maximum apparent velocity amplitude (defined to be \/2 times the rms value), 
and V is the value obtained by the excitation theory, scaled to render V = Vq for the Sun. The solar value was taken to 
be V0= 0.23 ms^^ (from Houdek & Gough 2002). 
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3.31 


60.00 
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18.13 


33.64 


9.04 
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"from Bouchy & CaiTier (2001), ''from Bedding et al. (2001), 
'^from Martic et al. (1999), ''Frandsen et al. (2002). 
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Figure 11. Predicted apparent velocity amplitudes (de- 
fined to be \/2 times the rms value) as function of light-to- 
mass ratio for stochastically excited oscillations in other 
stars. Observations from several authors are plotted by 
the plus and triangle symbols. The theoretical estimates 
by Houdek & Gough (2002) are plotted as diamond sym- 
bols. The scaling law by Kjeldsen & Bedding (1995) 
is illustrated by the dotted line and results reported by 
Samadi et al. (2005) are indicated by the dashed line. 

is the scaling law by Kjeldsen & Bedding (1995), and 
the dashed line is the scaling relation reported by Samadi 
et al. (2005) using the convective velocity profiles 
from numerical simulations (Stein & Nordlund 2001), a 
Lorentzian frequency factor in equation ( fT9t , and the the- 
oretical damping rates from Houdek et al. (1999). For 
hotter stars they find better agreement with observations. 

8. COMPARISON WITH NUMERICAL SIMULA- 
TIONS 

Recently Stein et al. (2004) reported stochastic excita- 
tion rates for radial p modes in various stars obtained 
from hydrodynamical simulations of their surface con- 
vection zones. By comparing the excitation rates between 
these simulation results and the model computations dis- 
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Figure 12. Acoustic excitation rates Pf of radial modes 
for the Sun and for the two solar-type stars a Cen A and 
Procyon. The solid curves are the results from numerical 
simulations by Stein et al. (2004) and the dashed curves 
are the model results of Sections 3 and 6. The model 
results have been scaled with a factor A such as to obtain 
the same maximum solar value of P{ than that from the 
numerical simulations (adapted from Stein et al. 2004). 



cussed in the previous sections, we can put some con- 
straints on the modelled oscillation amplitudes. In Fig. [12] 
the excitation rates for the Sun and for the two solar-like 
pulsators a Cen A and Procyon are plotted versus fre- 
quency. The solid curves are the results from the numer- 
ical simulations and the dotted curves are the excitation 
rates obtained from equation (fTsl l. The theoretical values 
were multiplied with a scaling factor A such as to obtain 
for the solar model the same maximum value of Pf than 
that from the numerical simulations. This allows us to 
compare the results for Procyon and a Cen A between 
the simulation (Psim) and model (Pmodci) computations 
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Figure 13. Total excitation rate Pf^tot and individual contributions from the turbulent Pf^turb and gas Pf,gas pressure 
fluctuations obtained from numerical simulations for a solar model (left panel) and for a model of an F star ( right panel ). 
For the solar model the two individual excitation sources are roughly similar in magnitude but for the hotter F star the 
driving from the Reynolds stress (turbulent pressure) fluctuations dominates. Partial cancellation between the Reynolds 
stress and gas pressure (buoyancy) fluctuations leads to a total excitation rate Pf _tot that is smaller (by a factor of about 
two in the F star, indicated by the two horizontal, dotted lines) than the contribution from the Reynolds stresses Pf,turb 
(adopted from Stein et al. 2004). 



(from here on we omit the subscript f). From this com- 
parison (see Fig.[T2]i we obtain for Procyon 

Pmodcl — 2 X Psim t (22) 

and from Table 2 we obtain for the ratio between the mod- 
elled (Kiiodoi) and observed iVohs) velocity amplitude for 
Procyon 

Kiodoi ^ 2.5 X Kbs . (23) 

With the assumption that the numerical simulations pro- 
vide the "correct" maximum value for the excitation rate 
P = Psim, we can estimate from the expressions V = 
^Pjlf]! (see equation (|6]l), (|22] | and (l23l l the amount by 
which the modelled values of Pmodoi and of (2r//)^^ are 
at variance with the observations: 

Knodol^V2xl.77xl4bs, (24) 

suggesting that the modelled excitation rate Pmodci is too 
large by a factor of about two and the modelled value 
of \J1t]I too small by a factor of about 1.77. In Sec- 
tion 4 we concluded that the modelled eigenfunctions, 
and therefore mode inertia /, reproduce the observations 
reasonably well. Consequently we can argue that it is pre- 
dominantly the theoretical damping rate 77 in \J2r\l which 
contributes to the factor of 1.77, i.e. 77 is too small (or the 
modelled mode lifetime too large). 
A similar analysis can be done with a Cen A, for which 
we obtain from Fig [12] and Table 2: 

Kiodci ^ ^ X 1.04 X Kbs , (25) 

suggesting that for this star the modelled damping rates 
77 are in good agreement with the observations (see also 
lower panel of Fig|9l). 

It is interesting to note that the numerical simulations by 
Stein et al. (2004) show partial cancellation between 



the two excitation sources T(u) and Q(s') (see equa- 
tion ( fT2] i) arising from the fluctuating turbulent pressure 
(Reynolds stresses) and gas pressure (buoyancy force); 
cancellation between acoustic multipole sources was dis- 
cussed before by Goldreich & Kumar (1990) and Osaki 
(1990). This effect is illustrated in Fig. [T3] which shows 
the total excitation rate Pf ^tot and the individual contribu- 
tions Pf.turb and Pf,gas from the turbulent and gas pres- 
sure fluctuations obtained from numerical simulations of 
the surface convection in the Sun and in a hotter star of 
spectral type F, the latter having model properties that 
are similar to that of Procyon. It is particularly striking 
that for the F star the effect of this partial cancellation 
leads, for the most prominent modes, to a total excitation 
rate Pf ^tot which is, on average, by a factor of about two 
smaller than the excitation rate from only the turbulent 
pressure fluctuations Pf ,turb (illustrated by the two dotted 
horizontal lines). One is therefore tempted to argue that 
the overestimated values of the modelled excitation rate 
Pf in Procyon (which is an F5 star) could be partially at- 
tributed to having neglected the gas pressure fluctuations 
in equation ( fT2] i and more importantly its cancellation ef- 
fect with the turbulent pressure fluctuations. Including a 
formulation of this cancellation effect into the excitation 
model of Section 6.1 could bring the estimated velocity 
amplitudes, particularly for hotter stars, in better agree- 
ment with the observations. 
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